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Abstract 



The depletion interaction between two parallel repulsive walls confining a 
dilute solution of long and flexible polymer chains is studied by field-theoretic 
methods. Special attention is paid to self-avoidance between chain monomers 
relevant for polymers in a good solvent. Our direct approach avoids the 
mapping of the actual polymer chains on effective hard or soft spheres. We 
compare our results with recent Monte Carlo simulations [A. Milchev and K. 
Binder, Eur. Phys. J. B 3, 477 (1998)] and with experimental results for the 
depletion interaction between a spherical colloidal particle and a planar wall 
in a dilute solution of nonionic polymers [D. Rudhardt, C. Bechinger, and P. 
Leiderer, Phys. Rev. Lett. 81, 1330 (1998)]. 
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I. INTRODUCTION 



In polymer solutions the overlap of depletion zones for monomers due to repulsive confin- 
ing walls or mesoscopic particles dissolved in the solution induces an important and tunable 
effective interaction potential Jl]. For example, this depletion interaction explains success- 
fully phase diagrams of colloid-polymer mixtures 0-fl- Recent experimental techniques 
facilitate even the measurement of the depletion force between a wall and a single colloidal 
particle In the context of such solutions confined to thin films and porous materials 

the geometry of two parallel walls has been extensively studied as a paradigmatic case |]8Hl4 



For strongly overlapping chains as realized in a semidilute polymer solution, chain flex- 
ibility is taken into account within self-consistent field theory or within the framework of 
phenomenological scaling theory P^T7|. On the other hand, in a dilute polymer solution 



different chains do not overlap so that the behavior of the polymer solution is determined 
by the behavior of a single chain. To a certain extent and under certain circumstances, 
a single chain can be modeled by a random walk without self-avoidance (ideal chain). In 



three dimensions this situation is closely realized in a so-called ^-solvent |fLqj . If the solvent 
temperature is below the #-point (poor solvent) the polymer coils tend to collapse |19| , p0 . 
However, in the common case that the solvent temperature is above the #-point (good sol- 
vent) the excluded volume (EV) interaction between chain monomers becomes relevant so 
that the polymer coils are less compact than the corresponding ideal chains. The emphasis 
in this work is on the latter situation and we investigate the effect of the EV interaction 
on the depletion interaction between two parallel walls as compared to the case of confined 
ideal chains. 

By focusing on long flexible chains in a system of mesoscopic size we obtain mostly univer- 
sal results which are independent of microscopic details [p^ , ^T| -p6|. Due to the universality 
of the corresponding properties it is sufficient to choose a simple model for calculating these 



results. For our investigations we use an Edwards-type model [[]J|,|n],^2j for the polymer 
chain which allows for an expansion in terms of the EV interaction and which is amenable 
to a field-theoretical treatment via the polymer magnet analogy. The basic elements in this 
expansion are partition functions zlcg(r, r') for chain segments which have no EV interaction 
(as indicated by the superscript [0]) and with the two ends of the segment fixed at r and r'. 
This perturbative treatment has to be carried out in presence of confining geometries. We 
consider two structureless parallel walls in d dimensions and a distance D apart so that in 
coordinates r = (rii, z) the surface of the bottom wall is located at z = and rn comprises 
the d — 1 components of r parallel to the walls. The surface of the upper wall is located at 
z = D. The interaction of the polymer with the non- adsorbing walls is implemented by the 
boundary condition that the segment partition function and thus the partition function for 
the whole chain vanishes as any segment approaches the surface of the walls |T8|j23 



i.e., 



4J(r,r')->0, z,z'^0,D. (1.1) 

For the present purpose the only relevant property which characterizes one of the polymer 
chains is its mean square end-to-end distance in the bulk solution which we denote by d 1Z 2 X 
|24^p6|; for convenience we include the spatial dimension d as a prefactor. The results 
presented in the following are obtained for d = 3 both for ideal chains and for chains with 
EV interaction. In Sec. [TJ we present our results for the interaction potential and the force 
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between two parallel walls. We note that for chains with EV interaction these results are 
only valid in the limit D S> TZ X because our theoretical approach is not capable to describe 
the dimensional crossover to the behavior of a quasi d — 1 dimensional system which arises 
for D <C TZ X . In Sec.|T| we compare these results with the simulation data of Milchev and 
Binder [H|. In Sec. |IV] we apply the Derjaguin approximation in order to obtain from the 



results in Sec. || the depletion interaction between a spherical particle and a wall in a dilute 
polymer solution and compare it with the corresponding experimental results of Rudhardt, 
Bechinger, and Leiderer ||. 

In view of the complexity of the actual experimental systems involving spatially confined 
colloidal suspensions dissolved in a solution containing polymers, in the past the corre- 
sponding theoretical descriptions relied on suitable coarse-grained, effective models and on 
integrating out less relevant degrees of freedom. The gross features of these systems can be 
obtained by mapping the polymers onto effective hard spheres as pioneered by Asakura and 
Oosawa |[27|| . In a more refined description Louis et al. f2l||2l| derived effective interaction 
potentials between polymer coils such that they behave like soft spheres. This approach 
allows one to capture the crossover in structural properties to semidilute and dense polymer 
solutions. Based on such a model Louis et al. |29| have calculated, inter alia, the correspond- 
ing depletion energy between two parallel repulsive plates. Besides presenting the depletion 
energy for ideal chains in terms of an expansion introduced by Asakura and Oosawa, they 
find the occurrence of repulsive depletion forces in an intermediate regime of D upon signifi- 
cantly increasing the polymer density. In our present completely analytic study we focus on 
dilute polymer solutions, for which the depletion forces turn out to be always attractive, by 
fully taking into account the flexibility of the polymer chains and the self-avoidance of the 
polymer segments with a particular emphasis on long chains. This complementary point of 
view allows us to make contact to the Monte Carlo simulation data in Ref. [|nj] and with 
the experimental data in Ref. [§; both these comparisons have not been carried out before. 



II. EFFECTIVE INTERACTION BETWEEN PARALLEL WALLS 

A. Grand canonical ensemble 

In a dilute polymer solution the interaction between N different chains can be neglected 
so that the total free energy of the system is N times the free energy of a single chain. We 
consider the polymer solution within the slit to be in equilibrium contact with an equivalent 
polymer solution in a reservoir outside the slit so that there is exchange of polymer coils 
between the slit and the reservoir. The free energy of interaction between the walls in such 
a grand canonical ensemble is given by 



z V z 



(2.1) 



where Z\\(D) is the partition function of one polymer chain in a large volume V containing 
the walls and Z is the corresponding partition function of one polymer chain in the volume 
V without the walls. In the thermodynamic limit V — > M d one has |24|- 



3 



z^vz h 



(2.2) 



with Zf, = J Rd d d r'Zt,(r, r') and where Zj,(r, r') denotes the partition function of one polymer 
chain in the unbounded solution with its ends fixed at r and r'. Correspondingly Z\\(r, r') 
denotes the partition function of one polymer chain within the volume V containing the 
parallel walls and with its ends fixed at r and r'. The volume V = Vi + Vo can be divided 
into the volume Vi within the slit and the volume Vo outside the slit. Since the polymer 
chain cannot penetrate the walls, whose lateral extensions are large, Z\\(r, r') is nonzero only 
if both r and r' are in Vi or in Vo so that 



Z\\{D)= [ d d r [ d d r'Z\\{x,r') = / d d rZ Q {z) + / dfrZ^z) 

JV JV Jv JVj 

with Zoj(z) = J Vq d r'Z\\(r, r'). In the thermodynamic limit V — > 
\n(Z\\/Z) in Eq. ( [2.1| ) can be expanded using Eqs. ( [2.2| ) and ( |2.3| ), i.e., 



(2.3) 

l d the logarithm 



In 



Zn(D] 



In 1 



Z\\{D)-Z 



\Z\\(D)-VZ h 



(2.4) 



V 



where the ratio 11 z is of the order V//V which tends to zero because the slit width D 

is fixed. Since in the thermodynamic limit the leading contribution to the first integral of 
Eq. 



Z h 



1 
V 



d d r 
UVo 



Zq{z) 
Z h 



1 ] + I d d r 

V/ 



ZAz) 
Z h 



is independent of the slit width D we find for Eq. (|27T 

ZAz) 



5F = —n p kbT 



d a r 



Z h 



d d r 



Zi{z) 
Z h 



(2.5) 



D=oo . 



with the number density n p = N/V of the polymer chains in the bulk solution. The sec- 
ond integral in Eq. (|2.5| ) reduces to the sum of two half-space (-£^5") integrals which yield 



contributions proportional to the area A of the walls p5 

Z H s{z) 



d d r 



V/ 



z h 



1 



D=oo 



2 / d d r 
'hs 



Z h 



1 



-2A 



Act 

n p k B T 



(2.6) 



where we have introduced the surface tension Act between the polymer solution and the 
confining wall [compare Eqs. (1.7) and (2.41) in Ref. |25|]]. Note that the mean polymer 
density within the slit is determined by the bulk density n p , i.e., in the grand canonical 
ensemble the chemical potential fi of the polymer coils is fixed instead of the number Nj of 
polymer coils in the slit (see Subsec. piBp . 

According to Eqs. ( |2.5|) and ( |2.6|) the total grand canonical free energy Q of the polymer 
solution within the slit, 



Q = —n p ksT ADuj 



(2.7) 



with the dimensionless quantity 
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^u D d: m , (2.8) 

U JO £>h 



can be decomposed as 

n 



Df b + 2f s + 5f . (2.9a) 



ripkBTA 

On the rhs of Eq. (gjgg) appear the reduced bulk free energy per unit volume 

fb = -1 , (2.9b) 
the reduced surface free energy per unit area 

and the reduced free energy of interaction 

5F 

Sf = —rr^ ■ 2.9d 

n p k B TA 

B. Canonical ensemble 

If instead of the chemical potential /x the number Ni of polymer coils in the slit is used 
as independent variable, the total free energy F of the polymer solution within the slit in 
the canonical ensemble follows from Q as the Legendre transform 

F(N I ) = S}(jjL(N I ))+ f ji(N I )N I , (2.10) 

where Q is given by Eq. ( |2.7|) . The chemical potential \i is related to n p via |30 



/x = k B T\n{n p A d ) , (2.11) 

where A is the thermal de Broglie wave length of the particles, i.e., polymer coils. Equation 
( [2.7D implies for dilute solutions 

dQ(u) —n. u d£l — 

iV/ = — = = . (2.12) 

dfi k B T dn p k B T 

Thus F(Ni) is given by 

F(N I ) = -k B TN I + k B TN I \n(^-A d \ , (2.13) 



with uj from Eq. 
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C. Free energy of interaction and mean force 



We employ the polymer magnet analogy |18| , |21| -p5| in order to calculate the partition 



functions Zhs(z), Zj(z), and Z\, as needed in Eqs. (|2.5|) , ( |2.6| ), and ( |2.g| ) for a single chain 
with EV interaction. They are functions of the parameter uq which characterizes the strength 
of the EV interaction and Lq which determines the number of monomers of the chain such 
that Lq equals 1Z 2 = 1Z 2 X / 2 for an ideal chain in the bulk, i.e., for u = 0. The well- 



known arguments of the polymer magnet analogy |IT|[2"T| -|2"5[ imply for the present case the 
correspondence 

Z||(r,r';L , A«o) = Ao-Lo^OO* i(r')> |^ =q (2-14) 

between Zu(r, r') and the two-point correlation function ($ 1 (r)$ 1 (r')) in a 0(N) symmetric 
field theory for an TV-component order parameter field $ = (<&i, ... , <&w) in the restricted 
volume V/. In Eq. (|2.14|) the operator 



C to ^ Lo = ^- I dt e Loto (2.15) 
J c 

acting on the correlation function is an inverse Laplace transform with C a path in the 
complex t -plane to the right of all singularities of the integrand. The Laplace-conjugate 
to of L and the excluded volume strength w appear, respectively, as the 'temperature' 
parameter and as the prefactor of the (<fr 2 ) 2 -term in the Ginzburg-Landau Hamiltonian 



V 



H{<f>} = I d d r{ ^(V$) 2 + + t^($T \ (2.16) 



which provides the statistical weight exp(— 7^{<fr}) for the field theory. The requirement in 
Eq. (|L1D describing the repulsive character of the walls imposes the Dirichlet condition 

*(r) = 0, z = 0,D , (2.17) 

on both walls. This corresponds to the fixed point boundary condition of the so-called 
ordinary transition PT| , P2"| for the field theory. For the renormalization group improved 
perturbative investigations we use a dimensionally regularized continuum version of the field 



theory which we shall renormalize by minimal subtraction of poles in e = 4 — d ||33|| . The 
basic element of the perturbation expansion is the Gaussian two-point correlation function 
(or propagator) ($j(r) <& 3 -(r'))[ ] where the subscript [0] denotes u = [see Eq. (|Al|) in 
Appendix HI]. 

The loop expansion to first order in Uq and the renormalization of u in Eq. (|2.8|) is 
completely analogous to that outlined in Subsec.IIA of Ref. |[25|| . Some key results of this 
procedure relevant for the present case are given in Appendix Ej. The final result for f s and 
5f on the rhs of Eq. fl2.9a| ) are given by Eqs. ( |A^ ) and ( |A7| ) in Appendix 0, where Af = 
for the present polymer case and e — 1 in d — 3. 

Figure|l| shows the universal scaling function for the free energy of interaction 

®(y) = ^-Jf ( 2 - 18 ) 
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with 5f from Eq. ( |2.9d| ) and the corresponding scaling function for the force 

m = (2.i9) 

in terms of the scaling variable 

y = D/n x . (2.20) 

Figure [1] shows both the behavior for ideal chains and for chains with EV interaction. For 
chains with EV interaction the present theoretical approach is only amenable to capture 
the behavior for D/1Z X large, i.e., i/ > 1. As expected the depletion potential and the 
resulting force is weaker for chains with EV interaction than for ideal chains, because the 
EV interaction effectively reduces the depletion effect of the walls. The inset of Fig. |l| shows 
both and V for ideal chains. The narrow slit limits read 

e(y^0) = --^ + y (2.21) 

'7T 



and 

r(!/-0) = -l . (2.22) 



Deviations from the linear behavior in Eq. Q2.21 ) become visible only for y > |. In the 



opposite limit y — ► oo the leading behavior for ideal chains is given by 



e(y^oo) = -4J-\e-y 2 / 2 (2.23) 

V 7rjr 



and 



T(y -> oo) = -A J -- e~ y2/2 . (2.24) 
V ny 

The depletion potential in terms of the scaling function Q(y) is attractive. But whether the 
bulk contribution D fb has to be taken into account in addition depends on the 'experimental' 
setup. In the case that the force is measured between plates immersed in a container filled 
with the dilute polymer solution such that solvent and polymer coils can freely enter the slit 
from the reservoir, only the scaling functions G and T are relevant. This is also true for the 
particle wall geometry discussed in Sec. |V|. But if no exchange is allowed as for the Monte 
Carlo simulation discussed in Sec.lTTl] the force /C defined in Eq. (|3.2|) is needed. 



III. COMPARISON WITH MONTE CARLO SIMULATIONS 

Monte Carlo simulations of polymers are well established both for ideal chains and for 
chains with self-avoidance. In this section we compare our results with the simulation of a 
polymer chain between two repulsive walls by Milchev and Binder which corresponds 



to the case studied theoretically here. These authors use a bead spring model for the self- 
avoiding polymer chain with a short-ranged repulsive interaction between the beads. 
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In Refs. || and |14J (see in particular Fig. 4 in Ref. |14| ) it is stated that the total force 
K, between the two walls is repulsive and diverges in the narrow slit limit, i.e., 

DK ( DY l/v . . 

where v = | for ideal chains and v = 0.588 for chains with EV interaction and lZ g is the 



radius of gyration of the polymer chain in unbounded space |l4j]. [For the definition of 1Z 



see Eqs. ( |4.2[ ) and fl4.3| ) in Sec. [TV].] The qualitative difference to the scaling function T(y) 



presented for the force in Subsec. [II Q is explained by the following arguments, (a) The total 
force between the walls is repulsive due to the contribution from the cost in free energy 
caused by the loss of the total available space for the polymer chains within the slit, i.e., 
the bulk pressure — /& contributes to the total force, (b) The total force /C diverges due to 
the change from the grand canonical to the canonical ensemble (see Subsec. |Il Bp . In the 



simulation in Ref. |14j] the number Nj of polymers in the slit is given by one polymer in the 
slit volume. Therefore we obtain the total force K, by differentiating Eq. ( |2.13| ) with respect 
to the slit width D and by setting Nj = 1, i.e., 



where D uo = —Q/(n p ksT A) is given by the rhs of Eq. ( |2.9a|) in conjunction with Eqs. (|A6|) 



and ( |A7| ) in Appendix The lhs of Eq. ( |3.2| ) corresponds to the quantity Df in Ref. |L4] . 



Figure^ shows the comparison of the simulation data of Ref. |14| and the corresponding 
theoretical result from Eq. ( |3.2| ), the latter both for ideal chains and chains with EV inter- 
action. Note that the theoretical curve for chains with EV interaction is only valid for large 
D/TZg. The curve for chains with EV interaction is closer to the simulation data than the 
curve for ideal chains, in agreement with the fact that the polymer chain in the simulation 
is a self-avoiding one. One possible reason for the remaining deviation might be that the 
chain in the simulation is too short to be fully described by the present field-theoretical 
approach. The deviation of the Monte Carlo simulation data at large values of D/lZ g from 
the power- law behavior at small values of D /1Z g (dashed line) occurs because tends to 
1 for large D /1Z g , which is not captured by Eq. (|3.1|) . 



Figures 5 — 9 of Ref. |14|] show density profiles for the simulated chains with EV interac- 



tion, including a comparison with the analytical result for the profile of ideal chains in the 
slit (Fig. 9 in Ref. fl^]). We want to mention here that the field-theoretical treatment of the 



monomer density for chains with EV interaction requires a perturbation expansion involving 
integrals over G * G * G * G instead of G * G * G (see Appendix [A]) and thus in view of the 
technical challenges is beyond the scope of the present study. Alternative information about 
the monomer density profiles beyond the ideal behavior can be found in Ref. p4j in which a 
self-consistent mean-field approximation is used to obtain the monomer density profiles for 
a single polymer chain between two repulsive walls. 



IV. COMPARISON WITH EXPERIMENT 

Rudhardt, Bechinger, and Leiderer have measured the depletion interaction between 
a wall and a colloidal particle immersed in a dilute solution of nonionic polymer chains in a 
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good solvent by means of total internal reflection microscopy (TIRM). They monitored the 
fluctuations of the relative distance of the colloid particle from the wall induced by Brownian 
motion. From the resulting Boltzmann distribution one can infer the corresponding effective 
depletion potential between the repulsive wall and the particle. 

In order to compare these experimental data with our results we apply the Derjaguin 
approximation [|35|| . In the limit that the radius R of the spherical particle is much larger 
than both 1Z X and the distance a of closest approach surface-to-surface between the particle 
and the wall, the particle can be regarded as composed of a pile of fringes. Each fringe 
builds a fringe- like slit with distance D = a + r^/(2R), where r\\ is the radius of the fringe. 
Thus the interaction between the particle and the wall is given by 

•**<«) .j,,^ r dvve (±- + v 4) , (4.1) 



n p k B T x J \TZ X 2 

where v is a dimensionless variable and Q(y) is the scaling function for the free energy of 
interaction for the slit geometry [see Eq. Q2.18Q 1. 

In Ref. || the interaction potential $depi(o) is measured for nonionic polymer chains in 
a good solvent for polymer number densities n p = 0, 7.6, 10.2, 12.7, and 25.5 /im~ 3 . All 
these polymer number densities represent a dilute polymer solution so that $dc P i( a ) is a 
linear function of n p [see Eq. (|2.1|) 1. Therefore our results are applicable and Fig. |3] shows 
§dc P \/n p as a function of a. The crosses in Fig. | correspond to those values of a for which 
the above mentioned linear behavior $depi oc n p (not shown) is in good agreement with 
the experimental data. Deviations from the linear behavior $depi oc n p for small and large 
particle-wall distances a (not shown) can be explained by the experimental method TIRM 
used in Ref. ||: a higher interaction potential implies a lower probability to find the particle 
at the corresponding distance causing lower accuracy. It turns out that the total interaction 
potential as the sum of depletion potential, electrostatic repulsion, and gravity is highest 
for short and large distances, which are those where the linear relationship $depi oc n p is not 
confirmed experimentally. Figure^] also shows the corresponding theoretical predictions for 
the experimental data, both for ideal polymer chains and chains with EV interaction, i.e., 
chains in a good solvent as realized in the experiment. Note that all parameters entering 
these theoretical predictions are fixed by available experimental data for n p , a, and 1Z g , i.e., 
there are no freely adjustable parameters. In particular, the radius of gyration 1Z g of the 
polymer chains used in Ref. has been measured fairly accurately by means of small angle 



scattering of x-rays 36], resulting in 1Z„ = 0.101/mi. According to the definition of 1Z„ as 



measured in small angle scattering experiments [p7H , i.e. 



T?2 _ Jd 3 rjd 3 r'p(r)p(r')\r-r'\ 2 



9 2 J d 3 r J d 3 r' p(r)p(r' j 

where p(r) is the monomer density, in d = 3 one has |38| 



7^ = 0.692772*, d = 3. (4.3) 

Figure^ shows that the experimental data of Ref. deviate from the theoretical result 
derived here. Note that this deviation is not visible if the polymer size is used as a freely 
adjustable parameter in order to gain agreement between the experimental and theoretical 
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data. The theoretical curve for chains with EV interaction, as realized in the experiment, 
is even further away from the experimental data than the theoretical curve for ideal chains. 
The latter observation confirms the necessity of reinterpreting of the experimental data, 
e.g., by adjusting the radius of gyration. This can be done such that agreement with the 
theoretical data for chains with EV interaction is gained for the intermediate values of the 
distance where the linear relationship $depi oc n p allows for a direct comparison with the 
theory. [The remaining differences outside this intermediate regime (o) pose a separate issue 
as discussed in the paragraph following Eq. (|4.1|) .] However the difference between the two 
theoretical curves for ideal chains and chains with EV interaction is small compared to the 
deviation from the experimental data. 



Any attempt to use the Monte Carlo data obtained in Ref. |L4| for predicting the deple- 
tion interaction between a colloidal particle and a wall would require two integrations of /C 
and a change to the grand canonical ensemble, which poses prohibitive accuracy problems. 
Nonetheless a qualitative statement can be given easily. Figure shows that the force ob- 
tained in the Monte Carlo simulation is weaker than the force calculated for chains with EV 
interaction. This is also expected to hold for the depletion interaction between a particle 
and the wall. 



V. CONCLUDING REMARKS AND SUMMARY 

Based on field-theoretical techniques we have determined the effective depletion inter- 
action between two non-adsorbing walls confining a dilute solution of long flexible polymer 
chains. Our main results are: 

1. The field-theoretical calculation yields the universal scaling functions of the depletion 
interaction potential and the corresponding force for ideal chains and for chains with 
excluded volume interaction in the limit y = D/1Z X ^> 1 , where D is the separa- 
tion between the walls and TZ X is the projected end-to-end distance of the chains [see 
Eqs. ( [2.i8| ), (|2.19|) , and (fO|) and Fig.[l|]. The depletion potential is weaker for chains 
in a good solvent than for ideal chains. 

2. For y > 1 we find fair agreement with corresponding Monte Carlo simulation data 
[H|j if the excluded volume interaction is taken into account (see Fig.^j). We surmise 



that remaining discrepancies are due to higher order terms in the field-theoretical 
calculation which are not yet taken into account and due to the possibility that the 
length of the simulated polymer chain has not yet reached the scaling limit for which 
the field theory is valid. 

Using the Derjaguin approximation we have compared our theoretical results with the 
experimental data P for the depletion potential between a spherical colloidal particle 
and a wall (see Fig. We obtain a fair agreement only if the radius of gyration TZ g of 
the polymers is used as a fit parameter. This value, however, differs from independently 
determined values for !Z g . The reasons for these differences are not yet understood. 
The excluded volume interaction between the monomers of the polymer chain plays 
only a minor role for reaching agreement between theory and experiment. 
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APPENDIX A: 



The Gaussian two-point correlation function in the slit of width D reads 
($,(!■) ^(rOho] = 8 ij G(r,r';t ,D) = 8 ij G(\r ll -r\ l \,z,z , ;t ,D) 



(Al) 



d^p 



(2 



7T 



— exp[ip • (r|| - rj|)] G(p, z, z; t , D) 



with the Gaussian propagator G(p, z, z'\ t , D) in p-z representation given by |32| , |39 

1 



G(p, z,z';t ,D) 



2b 



e -b\z-z'\ _ e -b(z+z>) 



+ 



,-b(z-z') _|_ e ~b(z'-z) 



-b(z+z') _ Mz+z') 



?2bD 



where b = \Jp 2 + t . The loop expansion of the total susceptibility reads |40 



D 



x(to,D,u )= / dzdz'G(0, z, z') 
Jo 

~Y~~T~ J dzdz ' dz " j (2 n y P i G( Q , z i z ") G(p, z", z") (5(0, z", z ) 
+ 0(ul) . 



(A2) 

(A3) 
(A4) 



The procedure outlined in Subsec. II A of Ref . [p^j yields the renormalized total susceptibility 
in one loop order as (see also Appendix B of Ref. [[51]) 

Xlw (T=(Dfj,)H,D,u) 1 2 1-e-v^ 



D 3 



t r 3 / 2 1 + 
M + 2 2 



(A5) 



3 r 3 / 2 



v^ + 2^- 



1 + e- 



-3 



1 - e^i 
1 + ev 77 



2/ x + 2 \n(Dfi) - In r + ln(47r) + 1 - C E - 8 / ds-^ 

i e 2 v TS — 1. 



-47T 



1-73 + e~^(2 - V3) + e- 2 ^(I _ ^ 



2 y/3 J 



1 + e 



<is- 



2 1 

+ 



1 + e-v 77 y o e 2 ^ s - 1 V s + i s - i ' s - 1 s + 1 



+ 0(« 2 
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Here \i is the inverse length scale which determines the renormalization group flow; t and u 
are the renormalized and dimensionless counterparts of t and u , respectively; Ce is Euler's 
constant and for the definition of the constant fi we refer to Ref. |25[ . 

The free energy is obtained via the inverse Laplace transform of Xren{D) and normal- 
ization by the transformed renormalized total susceptibility for the unbounded space Xren,&- 
The decomposition into bulk, surface, and finite-size contributions is carried out by the 
analysis of the scaling behavior of these parts of the free energy. The surface and finite-size 
parts of the free energy [see Eq. (|2.9| )1 at the fixed point of the renormalization group and 
for M = with the scaling variable y = D / 1Z X are given by 



1 - 



31n2 7T 



7T 

71 



(A6) 



and 



^ = 4Erf c m-4,/^e- 2 / 2 (A7) 

■ (4 + i| - 4vr + 6 \n(2y 2 ) - QC E ^j + Erfc(^=) (2n-^=-2 \n(2y 2 ) + 2C E ^j 
e^^lnr"! _ _ re _v/ ^lnr" 



-y 2 /2 



_yV2iT 



2y^ 



3£, 



-3/2 



+ 0(e 2 ) 



>BH\ 



where Xren in Eq. ( |A5| ) is expanded for large plate separations D up to order 0(e~ 
because for small distances the correct behavior including the dimensional crossover cannot 
be obtained even for the full expression. But using the expansion has the additional benefit 
to yield partly analytical results for the separations D of interest here. For d = 3, R x is 
related to the radius of gyration lZ g by Eq. QL3|) . The full result for ideal chains for 5f is 
given by Eq. for u = 0, i.e., 



5f 



-DC, 



1 

'2V 2 



1 



r 3/2 1 _|_ ev /7 



(A8) 



which is valid for all y. 
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FIGURES 
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FIG. 1. Scaling functions (a) Q(y) for the depletion interaction potential and (b) T(y) = — ^ 
for the depletion force between two repulsive, parallel plates at a distance D confining a dilute poly- 
mer solution in terms of the scaling variable y — D Jl^x [see Eqs. (|1JD and ( gig) ]. K x = V^TZg 
for ideal chains and TZ X = 1.444 TZ g in a good solvent for a bulk solution [see Eqs. ( |4.2| ) and ( |1.3| )] 
where lZ g is the radius of gyration. For ideal chains the full expression, which is valid in the whole 
range of y (solid line; see also the inset), and an expansion of this expression, which is valid for 
y > 1 (dashed line) are shown. The same expansion is shown for chains in a good solvent (dotted 
line). The dotted line stops where the dashed line starts to deviate appreciably from the solid line. 
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FIG. 2. Depletion force K, [Eq. (|3.2|)] between two parallel walls at distance D confining 
repulsively a dilute polymer solution. !Z g is the radius of gyration of the chains [see Eqs. ( |4.2[ ) and 
(|4.3| )1- The solid circles correspond to the Monte Carlo simulation data in Ref. [14]. The force 
obtained from Eq. ( p.2j ) is shown for ideal chains (solid line) and self-avoiding chains (dotted line). 
The latter line stops where the expansion for large D/lZ g turns out to become unreliable. The 
dashed line represents the asymptotic behavior at small distances for chains in a good solvent [see 
Eq. (3.1)], with a fit for the amplitude of the power-law behavior. For D/lZ g — > oo the reduced 
force DK,/ (ksT) tends to 1. 
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FIG. 3. Depletion interaction potential <t?depi( a ) between a spherical colloidal particle im- 
mersed in a dilute polymer solution and the container wall as a function of the distance a of closest 



approach surface-to-surface between the sphere and the wall [see Eq. (4.1)1. The interaction po- 



tential is given in units of fc^T and of the polymer number density n p . The circles and crosses 
indicate the experimental data from Ref. [6]. Crosses show the range where the linear relationship 
^depi oc %> allows for a direct comparison with the theoretical results. The theoretically calculated 
depletion interaction is shown for ideal chains (solid line) and chains in a good solvent (dotted 
line) as realized in the experiment. These curves correspond to the value TZ g = 0.101/mi, which 
has been determined by independent experiments [see Eqs. and (fD|)]. Using the radius of 
gyration as a fit parameter yields the dashed line corresponding to !Z g = 0.13/iin and ideal chains. 
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